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The SU{2) gluon propagator in Landau gauge is studied on the lattice. Our gauge fixing procedure 
employs simulated annealing and Z(2)-flips. It finds higher maxima of the gauge functional 
compared with those obtained with the standard overrelaxation and leads to systematic deviations 
of the gluon propagator in the infrared region. In particular, finite-size effects for lattice sizes 
from (1 .7 fm)"*^ up to (6.5 fm)^ become weak. The propagator shows a plateau at /? sa 300 MeV. 
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1. Introduction 

Over the years, considerable progress has been made in solving non-perturbative Dyson- 
Schwinger equations (DSE) for gauge-variant Green functions, in particular for the covariant Lan- 
dau gauge (for a recent review see Besides the interest in the DSE solutions as input for 
Bethe-Salpeter or Faddeev bound state equations, their infrared asymptotics is of importance for a 
check for gluon and quark confinement scenarios proposed by Gribov [Q] and Zwanziger [||] on one 
hand and Kugo-Ojima [Q] on the other. These scenarios claim confinement to be intimately con- 
nected with a Landau gauge ghost propagator diverging and with a gluon propagator vanishing in 
the zero-momentum limit. Such a behavior has been realized with asymptotic power-type solutions 
of (truncated) DSE with infrared exponents leading necessarily to a running coupling constant with 
a non-trivial infrared fixed point This behavior has been confirmed independently by studies 
of exact renormalization group equations ^ and with stochastic quantization [^. Recently it has 
been even argued that a unique and exact power-like infrared asymptotic behavior of all Green 
functions can be derived without truncating the hierarchy of DSE ||8]|. However, in order to inter- 
polate the full momentum dependence from the infrared to the perturbative ultraviolet regime, one 
still has to rely on truncations which are hard to control. Very recently, solutions of the truncated 
system studied on a finite torus have been presented with a specific finite-size dependence which 
smoothly turns into the exact power-like infrared behavior at infinite volume 

This gives us a good motivation to compare with the ab-initio non-perturbative path integral 
approach approximated on a Euclidean four-dimensional lattice. The lattice approach has its own 
limitations. Numerical simulations can be carried out only on a finite lattice. Therefore, for large 
momenta close to the inverse lattice spacing we shall encounter discretization effects, whereas 
at low momenta we are faced with the limitations of the finite volume as well as with rotational 
symmetry violations due to the hypercubic lattice geometry. Moreover, gauge fixing is not unique 
resulting in the so-called Gribov problem. It has been argued that the gauge copy dependence 
should disappear in the infinite-volume limit if the copies are bounded to the Gribov region - the 



positivity region of the Faddeev-Popov operator [jlOp. But on a finite lattice, Gribov copy effects 



may influence the infrared asymptotics and therefore, at least partly, be responsible for finite-size 
effects. Standard algorithms like overrelaxation (OR) find always local extrema of the gauge func- 
tional. Repeating such an algorithm with random initial gauges one can find better extrema coming 
closer to or eventually finding the global maximum - i.e. elements of the. fundamental modular re- 
gion. We call the copy found after a number of trials which guarantees stable values of the gauge 
functional, at least in the statistical average, the best copy (be) to be compared e.g. with the. first 
copy (f c). The gauge transformations determined in these gauge fixing procedures are normally 
restricted to be periodic. Under these circumstances for the Landau gauge the be ghost propagator 
has been shown to deviate up to 10% from f c results, whereas the gluon propagator did not differ 



within statistical errors [O, 12] 



In this contribution we present an improved gauge fixing method which allows to reach con- 
siderably higher extrema of the gauge functional than the above mentioned be OR method. In the 
simpler case of SU (2) we shall demonstrate the gluon propagator to become influenced and to have 
a weaker volume dependence in the infrared. We hope that this method will allow to check in the 
near future whether the gluon propagator really has the chance to tend to zero in the infrared limit. 
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The new method reUes first on a systematic use of the simulated annealing (SA) algorithm with 
subsequent finalizing OR to maximize the gauge functional and second on an enlargement of the 
gauge orbits by special non-periodic (modulo elements of the center Z(2)) gauge transformations, 
representing an exact (unbroken) symmetry of the local gauge action. First results obtained with 



this method were reported in [ |13[ |14[ ]. Restricting ourselves to the infrared region we present gluon 
propagator results here only for one considerably strong bare coupling value j8 = = 2.20. The 
corresponding lattice scale a is fixed with the string tension a = (440 MeV)^ adopting ^/aa = 0.469 
Thus, our largest lattice size 32"* corresponds to a volume (6.5 fm) . 

2. Improved gauge fixing 

Landau gauge fixing is equivalent to maximizing the gauge functional of a given lattice field {U} 

Fu[g] = ^fu[g], fu[g]=Y^\TrW,^ with 'U,^=g{x + fi)U,^g{xy (2.1) 

with respect to the local gauge transformation g{x) G SU{2). The SA method generates g 
stochastically with the Boltzmann weight w{g) oc exp{—fij[g]/T) , where the "temperature" 
T G [Tmin^Tmax] is a technical parameter which has to be lowered (we have chosen equal tem- 
perature steps between the lattice sweeps) from a certain value T,„ax until g is locked within the 
region of attraction of a local maximum. For the local updates of g the heatbath algorithm is used. 
After having reached Tmin, OR sweeps are employed until the lattice equivalent of d^Aj^{x) = 
is reached at all x with a given accuracy. The more slowly the SA cooling process is chosen the 
higher should be the probability to reach the global maximum. The method has been very suc- 
cessfully applied for the first time for gauge fixing in the case of the maximally AbeUan gauge in 



Ref. [|16|]. In order to see in as far the SA method (with finalizing OR) is more efficient than the 
only application of the OR algorithm we have selected for each gauge field up to 15 highest 

distinct local maxima F,, / = 1 , 2, On a lattice 16"* at j8 = 2.40 they can be well identified with a 

sufficiently large number of repetitions with initial random gauges. We measured the probabilities 
P{F) for each method to find the values F = Fi. The result is shown on the left of Fig. [l[ We 
compare OR with SA, the latter for various choices of the number of temperature iteration steps. 
SA is clearly seen to win even with a number of 0(1000) iterations. An SA iteration costs more 
CPU time than the simpler OR sweep. Therefore, one could think a repeated application of OR to 
be more time efficient. On the right hand side of Fig. |I| the probability to find the overall maximum 
Fi = Fffiax is shown versus the average CPU time required for the given version of algorithm. We 
see that even with respect to the computing time SA (including finalizing OR) is more efficient to 
find Ffnax- We are convinced that SA becomes increasingly efficient for lower jS and larger volume, 
respectively. Moreover, we have seen that SA is much improving, when microcanonical steps (in 
the following always three) are included after each iteration. In the left part of Fig. ^ we compare 
the probabilities P{F) of OR with SA for various lower temperatures r,„;„. We have measured the 



performance parameter introduced in Ref. [17] and defined as G = — log(l — P{Fmax))/tr, where 
tr denotes the CPU time in arbitrary units. Corresponding estimates for G are shown in the right 
part of Fig. ^. For what follows we decided to apply SA with r„„>, = 0.01 and 1000 iterations with 
equal temperature steps always combined with microcanonical sweeps. 
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Figure 1: Left: Probability to find the gauge copy with the functional value Fi of rank / = 1,2, .. . for 
the SA method with a number of temperature steps varying from 1000 to 6000 (for fixed T,ni„ = 0.4 and 
Tmax = 1-4). For comparison we show also the result of the standard OR method with one (the "first") 
random copy. Right: The corresponding probability to find the overall highest maximum F\ = F,„ax is 
shown vs. CPU time required for the SA method with varying number of temperature steps. For comparison 
the result for the OR method is shown when repeated with an increasing number of initial random gauges 
(curved line). The CPU time unit is the average time the OR method needs for one gauge copy to achieve 
the required accuracy of gauge fixing. 37 configurations, each with 50 gauge copies have been considered. 



The second feature of our improved gauge fixing procedure are Z(2) flip transformations. 
For SU (2) gauge theory, each flip transformation consists of a simultaneous Z(2) flip of all links 
Ux^ —>■ — Ux^ throughout a 3D hyperplane at a given value of the coordinate x^. This is just a 
particular case of a gauge transformation which is periodic modulo Z(2), 

g{x + Lfi)=z^g{x), z^ = ±lGZ(2). (2.2) 

The procedure is equivalent to search for the best sector (determined by the signs of the four aver- 
aged Polyakov loops) among 2^ = 16 sectors that provides the highest maximum of F. In order 
to decide which one is the optimal sector, the SA method has to be applied repeatedly with the aim 
to find the best copy within each sector. In practice the procedure can be somewhat simplified (see 



Ref. []14j]). We will abbreviate the combined gauge fixing method as the FSA (flip-SA) algorithm. 

In Table[T]we show the strong effect of the flips on the average gauge functional. F{nc) denotes 
the best functional value found with SA from ric random starts in every chosen flip sector. In case 
"SA" we did not apply flips at all, i.e. the Polyakov loop sector is chosen randomly by the Monte 
Carlo procedure. In case "FSA" we have searched within all 16 flip sectors. 



3. Lattice gluon propagator results 

The combined FSA method, and for comparison also the standard OR method, have been applied to 
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Figure 2: Left: Probability to find the gauge copy with the functional value Fj of rank ; = 1,2, . . . , 15 
shown for the methods OR and SA, for the latter with varying final temperature r,„,„ — 0.01,0.2,0.4 and 
a fixed number of 1000 temperature steps each supplied with 3 microcanonical steps {T,nax = 1-4). Right: 
Performance parameter G as defined in the text for the S A method shown as function of the final temperature 
Tmin- For both figures 33 configurations, each with 50 gauge copies have been considered for j3 = 2.4 and 
lattice size 16'*. 
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Table 1: The average gauge functionals {F{nc)) with an arbitrary value = 0.82800 subtracted. For 
the lattice sizes 16^ and 24^ the numbers of investigated MC configurations with j3 — 2.20 are 60 and 46, 
respectively. 



the computation of the gluon propagator at momentum = {2/a) sin {nk^/L), G (— L/2,L/2] 
flb _ /7« tiAl^bt ,A\ _ / X '^''\5"^D{p) , (3.1) 



(3.2) 



D"^%{p) = {A^^{k)A>',{-k))=[5,, 



where A{k) represents the Fourier transform of the gauge potentials 

1 



A^(x + /x/2) 



after the gauge has been fixed. 
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Figure 3: Gluon propagator versus momentum and the zero-momentum propagator D{Q), for /3 = 2.20 and 
for various lattice sizes, obtained with f c OR compared with be FSA. The lattice size is 32^. 
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Figure 4: The gluon propagator and £>(0), obtained with be FSA gauge fixing, shown in the infrared region 
for various lattice sizes (j3 = 2.20). 



Fig. m shows the comparison of the f c OR results obtained for several lattice sizes with the 
be FSA result for 32^ only. At = the zero-momentum data points D(0) are also plotted. The OR 
data exhibit quite strong finite-size effects. Contrary to the OR results the FSA data seem smoothly 
to extrapolate to the D{Q) data point. In Fig. ^ we show our be FSA result for various lattice sizes. 
In comparison to OR (see Fig.^ the FSA result shows considerably less finite-size effects down to 
the lowest accessible momenta. All data points fall more or less onto a universal curve. This leads 
us to hope that the visible plateau indicates the existence of a turning point beyond which D{p) 
starts to decrease for p — > 0. 
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4. Conclusions 



In this contribution we have discussed an improved gauge fixing method which takes Z(2) flips 
into account and makes consequently use of simulated annealing to maximize the Landau gauge 
functional. The combined algorithm finds considerably larger functional values. It lowers the 
values of the gluon propagator in the infrared in comparison with the OR results. Moreover, finite- 
size effects seem to become suppressed. They do not show the specific behavior found with DSE 
on a finite torus [^]. By further increasing the lattice size we hope to see D{p) to pass a maximum 
and to tend to smaller values in the far infrared. So far, such a behavior has been found only in the 



lower dimensional cases [18, 19]. 
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